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Abstract. We show that, above the critical temperature, if the dimension D of 
a given Ising spin glass model is sufficiently high, its free energy can be effectively 
expressed through the free energy of a related Ising model. When, in a large sense, 
D = oo, in the paramagnetic phase and on its boundary the mapping is exact. 
In this limit the method provides a general and simple rule to obtain exactly the 
upper phase boundaries. We provide even simple effective rules to find crossover 
surfaces and correlation functions. We apply the mapping to several spin glass 
models. 
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1. Introduction 

The spin glass model represents one of the most difficult challenges in physics (see 
Ref. m and referred articles). When compared to a non random model, the spin glass 
model turns out to be hugely more difficult to solve. In fact, to solve a spin glass 
model, one should be able to study, at least in principle, either a non random model, 
but with non uniform frozen couplings, or, as appears evident in the framework of the 
replica approach, a system of n coupled uniform non random models, extrapolating 
some suitable derivative in the limit n — )■ 0. From this point of view it seems fairly 
hard that a spin glass might be solved simply in terms of a mapping with a non 
random and uniform model. Nevertheless, in this paper we show that, whereas in the 
most rich and complex low temperature phases such a mapping is impossible, when 
the dimensionality D of the model tends to infinite, in the paramagnetic phase and 
on its boundary, the given spin glass model can be easily solved by simply considering 
a suitable corresponding non random Ising model: "the related Ising model" . More 
precisely, we shall show that, when in a large sense the dimensionality D is infinite, as 
happens, e.g., but not only, in mean field models and in generalized tree-like structures 
(allowing also for the presence of loops), once the phase boundary of the related Ising 
model is known, it can be immediately used in the spin glass model to find exactly 
the region of the paramagnetic phase (P) and its boundaries with the other phases: 
ferromagnetic disordered (-F), spin glass (SG), or antiferromagnetic disordered {AF). 

The derivation of this mapping, "spin glass" — > "related Ising model" , will be 
obtained in the framework of a replica approach, which however differs from the 
standard one j^. As we shall show, the use of replicas in the high temperature 
expansion of the free energy leads to a different procedure in which there is no 
functional to be extremized and, furthermore, when D is large, at high temperature 
a simple combinatorial approximation applies which, in the limit D —^ oo, becomes 
exact. We point out that, unlike the standard replica approach, in the proof for 
this mapping we do not need to rely on any ansatz concerning the choice of finding 
stationary points. Therefore our proof, up to question of the analytic continuation n 
integer — *■ n real, is exact 

Previous uses of the high temperature expansion to study spin glass models go 
back to [7] for the Sherrignton-Kirkpatrick model (SK) [S], and to [Sj and ^U] for the 
Ising spin glass models on hypercubic lattices with simmetric disorder, where very 
accurate results were found for small and high dimensions in [Hj and ^01 , respectively 
(we cite also the Ref. and the Refs. jElEl for the Potts glass model). However, 
as will be shown, our use of the high temperature expansion of the partition function 
is quite different. We do not use this expansion to directly calculate the averages over 
the disorder of a specific model, we use instead the expansion to find a general link 
between the spin glass and a suitable related Ising model. Once this link has been 
established, the singularities of the related Ising model will provide the singularities 
of the spin glass, signaling the phase transitions in the sense of Lee and Yang |14| . 

The practical potential applications of this mapping are remarkable. In fact, 
though the mapping may be used to analyze only the regions at high temperature, 
including the upper phase boundary and, within a certain approximation, the crossover 
surfaces and the correlation functions, as will be clear, its generality and simplicity 
make the approach particularly suitable to face those infinite dimensional spin glass 
models whose great complexity, e.g. due to the presence of too many parameters, 
can make difficult, or avoids at all, the application of the standard methods for spin 
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glasses even in the easier paramagnetic phase. In our approach instead, the model 
to be solved is the related Ising model which, as a non random model, turns out to 
be hugely easier to solve, analytically and/or numerically. In this work (part I) we 
will show this in several non trivial known examples, whereas in a forthcoming paper 
(part II) the mapping will be applied for studing the Ising spin glass model on general 
graphs and networks. 

The paper is organized as follows. In Sees. 2 and 3 we introduce the models 
and the condition under which the mapping holds. In Sec. 4, after introducing the 
definition of the related Ising model, we present formally the mapping, subdiving the 
results in the subsections 4.1 (same disorder for any bond), 4.2 (generalization), 4.3 
(upper line of the phase diagram). As a by-product, in the further subsection 4.4, 
via analytic continuation, we, improperly, force the mapping to find crossover surfaces 
and correlation functions. We stress that this extension of the mapping is not exact, 
but provide however a first effective insight about the physics of the model. 

The rest of the article until Sec. 11 is devoted to the proof of the mapping. In Sec. 
5 we give some preliminary ideas, whereas in Sec. 6 we recall the high temperature 
expansion for a general Ising (even non uniform) model. In Sec. 7, starting form the 
high temperature expansion, we carry out formally the average over the disorder to be 
used in Sec. 8 for calculating the free energy. In Sec. 9 wc specialize the expansion for 
a centered measure, whereas in the following subsection 9.1 the basic approximation 
in high dimension will be introduced and used to find the mapping for a hypercube 
lattice. In Sec. 10 and in the following subsection we generalize the mapping to any 
measure; finally, in Sec. 11 we extend the mapping to any model whose dimension, 
in a large sense, turns out to be infinite, as happens, in particular, but not only, in 
tree-like structures and generalized tree-like structures. 

The sections 12 and 13 arc devoted to some applications in the case of D finite 
and infinite, respectively. We anticipate that our approach takes into account only 
the leading term of a 1 / D expansion, so that the applications in the finite dimensional 
case are basically meant to show how the method works. 

Finally, conclusions and some outlooks are reported in Sec. 14. The paper is 
equipped with three appendix. 

2. Models 

Let us consider a D dimensional hypercube of side L, A = {1, . . . , L}^, and the set of 
the bonds b connecting two first neighbors sites 

T = {h= (ib, ib) : ib, jb e A, ib and jb first neighbors, ib < jb}- (1) 

We will indicate with N the total number of points, N = . More in general, we 
shall consider also systems over a graph. Given a graph g of N vertices, the set of 
links will be defined through the adjacency matrix of the graph, gi j = 0, 1: 

T = {h= {ib,ib) : ib,jb G .9, g^j, = 1, h < jb}- (2) 

The set of links of the fully connected graph will be indicated with Tf-. 

Tf = {b ^ {ib,ib) ■■ ib,jb ^ i, - ■ ■ ,N, ib < jb}- (3) 

The Hamiltonian of the spin glass with two-body interactions can be written as 

N 

H {{<J,}; {Jb}; {h,}) = - ^ Jbab ^ h,a,, (4) 
ber 1=1 
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where the hiS, are arbitrary external fields, the Jf,'s are quenched couplmgs, Ci is 
an Ising variable at the site i, and o-f, stays for the product of two Ising variables, 
= (^ib'^jbi with ill and jb such that b = {ib,jb)- 
The free energy F is defined by 



{~J^/J') , (8) 



- /3F = J dV {{Jb}) log {Z ({ JJ; {h,})) , (5) 
where Z {{Jb}', {hi}) is the partition function of the quenched system 

Z{{Jb};{K})^Y.^'^"^^'''^''^-''^''^''^^' (6) 

and dV {{Jb}) is a product measure over all the possible bonds b given in terms of 
normalized measures dfXb ^ (we are considering a general measure d/ib allowing also 
for a possible dependence on the bonds) 

dV {{Jb}) = n '^f'b (-^b) , / dfib {Jb) = 1. (7) 
beFf 

Among the measures most considered in literature, we cite in particular the 
Gaussian measure 

dfi{Jb) 1 

= , ^ exp 

""Jb V27rJ2 

and the "plus-minus" measure 

^-ljiJ-Jb) + ls{J+J,), (9) 

where and J represent disorder parameters. We will take the Boltzmann constant 
Kb = 1 . a generic inverse critical temperature of the spin glass model, if any, will be 
indicated with /3c; sometimes we will use the symbol < A > for the quenched thermal 
average of the quantity A (fixed values of the couplings {Jb})', finally the density free 
energy in the thermodynamic limit will be indicated with / = /(/?) 

/(/?) = hm F{(3)/N. (10) 
3. Dimensionality - Condition for the mapping 

Our mapping will be first derived in a D dimensional hypercubic lattice by applying 
an approximation becoming exact in the limit D oo. Successively, we will consider 
more general structures, whose dimension, in a large sense, goes to infinite as well. For 
a hypercube lattice the dimension D is related to the number of first neighbors of a 
vertex, 2D, so that, in this case, and in others in which the numbers of first neighbors 
is proportional to _D, _D oo if the number of first neighbors goes to infinite For these 
cases, we shall say that the set of links F is infinite dimensional in a strict sense. For 
other structures, a usual definition of dimensionality is 

;^oo log {I) 

where N{1) = 1 + mi + . . . + m; represents the total number of vertices within I steps 
of an arbitrary fixed root vertex 0, being the number of vertices at distance k 
from the root. We will not make use of this definition for the dimension, but here we 
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just observe what follows. When apphed to a hypercube lattice, Eq. l(TT]l returns the 
true dimension of the lattice D, whereas in the case of a tree with an average degree 
greater than 1, Eq. Ullfl gives D = oo. As an heuristic argument we can say that, 
given a structure F, the hypercube lattice most similar to F must have a dimension 
D given by Eq. Hll() : therefore, if, for F, D of Eq. goes to infinite, the mapping 
must hold. However, we will not need to consider this heuristic argument. In Sec. 
11, we will work out directly the derivation of the mapping for tree- like structures, 
little differences being involved with respect to the hypercubic lattice. Furthermore, 
we shall now introduce the most general condition under which the mapping is exact, 
stressing that the tree-like structure is not necessary. 

Let us define a path as a succession of first neighbor vertices connected by different 
bonds of F. A bond has length 1. The paths can be open or closed (loops). As will 
be clear later, depending on F and near the limit f3 ^ /3~ , the paths giving a finite 
contribution to the partition function Z arc infinitely long paths which can be: closed, 
as for hypercubic lattices; open, as for tree- like structures; and both closed and open, as 
for generalized tree-like structures. We say that F has a generalized tree-like structure 
if for any vertex there is at most a finite number of loops. 

What will emerge in subsection 9.1 and Sec. 11, is that 

when, in the thermodynamic limit, for the given set F, the total number of infinitely 

long paths per vertex goes to infinite and, choosing randomly two of them, the 
probability that they overlap each other for an infinite number of bonds goes to zero, 
in the limit (3 f3^ , the mapping becomes exact. 

We shall say that such a F is infinite dimensional in the large sense. It is immediate to 
verify that in the case of tree-like structures the above condition is trivially satisfied 
and similarly for generalized tree-like structures. In general, given F, the absence 
of loops, or the presence of a finite number of loops per vertex, turns out to be a 
sufficient condition for satisfying the above requirement. Note however that this is 
not a necessary condition. Loops can also be massively present; what is necessary 
is only that loops overlap, cither each other or with other paths, only partially (see 
Fig. Of course such a behavior does not happen in a hypercube lattice of finite 
dimension D. It is in fact easy to see that, in this case, the total number of paths 
of length I per vertex grows as c{l) ~ exp(/i£)/), with a suitable constant |21| . 
whereas the maximum number of paths which after I steps do not share any other 
bond grows only as Ci{l) ^ l^ . Therefore, if D is finite, in the limit I oo, we see 
that the condition cannot be satisfied. 

4. The mapping 

Given a spin glass model trough Eqs. H2I7|) . we define, on the same set of links F, its 
related Ising model trough the following Ising Hamiltonian 

TV 



where the Ising couplings 's have non random values such that V 5, 6' G F 




(12) 




or 




(13) 



(14) 
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Figure 1. Example to clarify the condition given in Sec. 3. Prom the root 
vertex depart several infinitely long paths. For simplicity, in the figure it is 
understood that, except for the chain indicated by an arrow, besides the loops 
a, b, c, d, L and L', there are no other loops connected to the root 0. It is also 
understood that some branching of the paths keeps on also out of the circle. From 
each of the loops a, b, c, and d, and from the other branches without loops, pass 
an infinite number of infinitely long paths, but the probability that two of them, 
randomly chosen, overlap for an infinite number of bonds, goes to 0. Through the 
chain of loops (infinite chain) pass also an infinite number of infinitely long paths, 
but in this case, the probability that, two of them randomly chosen, overlap for 
an infinite number of bonds is not 0. A way to generate, from the above graph, 
a graph infinite dimensional in the large sense, consists in deleting, e.g., all the 
infinite bonds of the chain indicated with dashed segments. The loops L and L' 
are two example of infinite closed paths. Note that, since they share only a finite 
number of bonds each other, their presence do not alter the condition of Sec. 3. 



In the following a suffix / over quantities such as Hj, Fi, //, etc. . . , or I3c^\ 
etc. . . , will be referred to the related Ising system with Hamiltonian 112(1 . We can 
always split the free energy of the quenched model as follows 

N 

-PF = Y, dfit, log (2 cosh(/3 Jf,)) + log (2 cosh(/3/i,)) + 0, (15) 
ber-^ i=i 

4> being the high temperature part of the free energy, and similarly for the related 
Ising model. Let (p be the density of (f> in the thermodynamic limit 

tp= lim (j)/N, (16) 

and similarly for ipi , the high temperature part of the free energy density of the related 
Ising model defined through Eqs. H12I14|I . As is known, the high temperature part (p 
(ipi) can be expressed in terms of the quantities Zf, = tanh(/3Jf,) (z^^' = tanh(/3J^'^'')) 
and Zi = tanh(/?/ii), i.e., the adimensional parameters of the Ising high temperature 
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expansion: 

if =<^({tanh(/3Jb)};{tanh(/3/i,)}), (17) 

ipi = ipi ({tanh(/3J<'^)}; {tanh(/3/i,)}) , (18) 

Hereafter, if not explicitly said, we set hi = 0. By varying in the Ising Hamiltonian 
l|12|l the couplings J^'^' with the constrains H13I14|I . we explore the function ipi and 
the non analytic points of ipi, if any, will signal some phase transition. The critical 
behavior of (pi will be characterized by an equation of the type Gi{{zl^^}) = 0, whose 

solution is provided only by universal, i.e., not depending on the {J^^^^}, quantities. 
A point on the critical surface S/, solution of the equation G/ = 0, will be indicated 
with {wlf^}. The surface E/ represents the boundary of some domain f^/ inside which 
the high temperature expansion providing tpi converges. As we shall discuss later, the 
domain ^/ is a convex set so that, if l^^^^l = | tanh(/3J^^^)| < for any 6, then 

Our main result concerns and its singularities, i.e. the phase boundaries of 
the spin glass model. We find convenient to separate these results in the cases of 
homogeneous and inhomogcneous measures. 

4.1. Case of an homogeneous measure (same disorder for any bond) 

Let be dfii, = dji for any bond 6 of F. so that, according to Eqs. 1)13(1 and 1)14(1 . 
the related Ising model corresponds to a homogeneous Ising model having a single 
coupling J^^^-* = and its critical behavior will be characterized by, at most, two 

points w^P > and < 0, if any. 

First, let us consider a system having a number of first neighbors per vertex 
proportional to D. Let be pi^'^^ and I3i^^^^\ respectively, the solutions of the 
equations, if any 

dfi tanh^ {[if ^'^Jb) =^wP, D>2, (19) 

d^i tanh(/3f /^^) Jb) = wP/^p, D>1, (20) 

where F or AF, in the l.h.s and r.h.s. of Eq. 1(20(1 . are in correspondence. We will 
show that, asymptotically, at high dimensions D, the critical inverse temperature of 
the spin glass model, /3c, is given by 

;9. = min{/3(^«),/3f/^^)}; (21) 

and in the paramagnetic phase the following mapping holds 

where, g{x) is some bounded continuous ftmction of order 1 such that g{x) for 
a; ^ 00; and p^eff, is given by 

P>eff = yV'/ (^j dii tanh'(/3Jb)^ , D > 2^-\ (23) 

where 

r 2, if 1(^7 (/d/x tanh2(/3Jb)|) >2|^7(/d// tanh(^Jb)) I, 

\ 1, if \ipi {J dfi tanh^ {(3 Jb)\) <2\(pi {J dfi t&nh{(3Jb))\. ^ ' 



= 0('3(r_^l V/3</3„ (22) 
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If the system F is infinite dimensional only in the large sense (Sec. 3), in the above 
equations D can be settled as infinite and, in the limit /3 ^ , Eq. holds exactly. 

4-2. Generalization 

The generalization of the above mapping to the case of an arbitrary (not homogeneous) 
measure dfib, useful for example for anisotropic models, in which wc have to consider 
even a given number of bond dependencies, follows straightforward. In this case, the 
related Ising model is defined by a set of, typically few, independent couplings {J^,^^}, 
trough Eqs. H13I14|I and its critical behavior will be fully characterized by the points 
of S/, solution of the equation G/ ~ 0. Equations (|19I24() arc generalized as follows. 

First, let us consider a system having a number of first neighbors per vertex 
proportional to D. Let be pi^'^'^ and I3c^^^^\ respectively, solutions of the two set of 
equations (if any) 

Gi ({ j dub tanh2(/3f<=) Jt)}^ =0, D > 2, (25) 

Gi (^{ Jd^lt tanh(/3f J,)}^ =0, D > 1. (26) 

Asymptotically, at sufficiently high dimensions D, the critical inverse temperature of 
the spin glass model f3c is given by 

/3. = min{/3P«',/3fM^)}; (27) 

and in the paramagnetic phase the following mapping holds 

-of^i^li^V V/3</3., (28) 



D 



Veff = jVi (^{ J dfib tanh'(/3J6)}^ , D > 2'-\ (29) 



where 

r 2, if \ipi ({/ d^ib tanh2(/?Jb)}) I > 2\ipi ({/ d^ib tanh(/3Jfc)}) |, , . 

\ 1, if ({/ d^lb tanh2(/3Jb)}) I < 2\^i ({/ dfit tanh(/3J,)}) |. ^ ' 

If F is infinite dimensional only in the large sense, in the above equations D can be 
settled as infinite and, in the limit /3 — > /3~ , Eq. holds exactly. 

4-3. Phase diagram: upper critical surface 

If D — > oo in the strict sense, Eqs. (|21I24|I and their generalizations, Eqs. H27I3()|I . give 
the exact free energy in the paramagnetic phase (P) ; the exact critical paramagnetic- 
spin glass (P - SG), and paramagnetic-P/AF (P - F/AF), (3f'^''\ surfaces, 
the stability of which depends on which of the two ones is the minimum. In the 
case of a homogeneous measure, the suffix F and AF stay for ferromagnetic and 
antiferromagnetic, respectively. In the general case, such a distinction is possible 
only in the positive and negative sectors of E/, whereas, for the other sectors, we 
use the symbol F/AF only to stress that the transition is not P — SG. Finally, 
the constrains > 1 or _D > 2 along the equations stress that the mapping is not 
otherwise defined. If instead D ^ oo only in the large sense. Eqs. H21I24|) and their 
generalizations, Eqs. H27I30|I . arc in general exact only in the limit (3 ~> (3~ and give 
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the exact critical paramagnetic-spin glass {P — SG), (3c , and paramagnetic-F/AF 
(P - F/AF), I3i^'^^\ surfaces. 

Notice that at zero field, due to the inequality 

^<Vi{{j dfib tanh(/3Jfc)}^ , (31) 

easily derived from the convexity of the logarithm function and, due to the mapping 
l|28|l . at least for centered measures, our estimations for the free energy become trivially 
zero if 13 — > oo in the strict sense: 

lim V3 = 0, for (3< (3c, (32) 

D^oo 

and the basic role of Eqs. 1)22124(1 . or ((28I3()() . is that to show how, in this limit, if 
approaches zero and which are its singularities. 

We stress that the presented result is obtained in a replica approach method but 
does not rely on any ansatz concerning the choice of finding stationary points. In 
known examples we have so far considered, Eqs. (|21I24() give results coinciding with 
those obtained in the framework of the standard replica approach equipped with the 
replica-symmetric ansatz, generally accepted as exact above the critical temperature. 
We recall that, above the critical temperature, the replica symmetric solution has 
been proved to be the maximum of the functional appearing in the standard replica 
approach, only in the SK model and in the "p-spin" models We stress also 

that the mapping is exact for any case in which T is infinite dimensional and with 
arbitrary measures including, for example, generalizations of the SK model, Bethe 
lattice models, and models defined on generalized tree-like structures. Finally, Eqs. 
H27I30|I are useful for more general models, like anisotropic models and models defined 
on bipartite lattices. 

Several general properties can be immediately derived from the critical conditions, 
Eqs. H19I23|1 . or their generalization. Eas. I(25I29|I . In particular we find as corollaries: 

i) A critical (3^^'^^ exists and is finite if f the related Ising model has a finite 
critical /3c^^ 

ii) If an Ising model on T with a homogeneous coupling J > has a phase 
transition at some /?c^^ , then the spin glass on T with a homogeneous "plus-minus" 
measure d^{Ji,)/dJb = S{J — Jb)/2 + S(J + Jb)/2, will have a SG transition for 
(3^'-^ solution of 

tanh2(/3(^^'j) =tanh(/3(^)j), (33) 
and, as a consequence, for the plus-minus measure one has always 

/^f^) > /3(^) (34) 

in) A critical /?c , exists and is finite if f for the related Ising model exists a 
point {w^^^} e S/ such that for wl^^ > one has 1 — 2 dfibiJb) > and 
for w^-^' < one has 1 - 2 dfibW < w[^\ ybeV 

iv) In the space of the parameters of the probability distribution d/i (for simplicity 
here we consider only the case of a same disorder for any bond) , possible lines of 
coexistence SG — F or SG — AF, must intersect the multicritical points satisfying 
the systems of equations given respectively by 

I J dn tanh^ {PcJb) = w''p\ ^^^^ 
1 J d^i tanh(/?cJfc) = Wp^; 
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/ d^i ianh^ {(icJb) = w^^^ 



£ ' (36) 



Analytic continuations: coexistence surfaces and correlation functions 

The statements until now presented are exact. In this subsection wc inchidc the 
following by-product. For D = oo, below the critical temperature, Eqs. (|22I24|) or 
H28I30|I are no more valid. This limit has the peculiar feature that, despite the fact 
that above the critical temperature the free energy is exact, in general, it cannot be 
analytically continued to lower temperatures ^Hl and this happens due to the fact 
that when D = oo, below the critical temperature the related Ising model is ill defined 
or, in other words, its thermodynamic limit does not exist; the density energy being 
infinite for any non zero value of the mean magnetization. 

What wc argue, instead, is that the analytic continuation of some physical 
quantities below even if it is not rigorous, provides a certain effective 

approximation. This is in particular the case for the coexistence equations H35t 
136(1 for evaluating the spin glass- ferromagnetic, {SG — F), and the spin glass- 
antiferromagnetic, {SG — AF), boundaries. As we shall show in the examples we 
will consider, these crossovers are, roughly speaking, close to the Almeida Thouless 
lines. Due the fact that the free energy provided by the mapping is eaxct in all the 
paramagnetic phase as Z? — *■ oo in the strict sense, we expect that the above analytic 
continuation turns out to be better in these cases, rather than in the cases in which 
D oo only in the large sense. 

Similar comments hold for extending the free energy part if to the case of arbitrary 
external fields hi. As we will discuss in Appendix A, the natural extension for ip 
consists simply in adding the furhter set of arguments {tanh(/3/ii)} as 

V^eff = jfi ({ J dfi tanh'(/3J6)}; {tanh{l3h,)}^ , (37) 

where I is to be determined with the analogous of Eqs. ((24|l . 

Now, if in infinite dimensions this extended mapping holds, due to the 
arbitrariness of the external fields hi in ip, by derivation wc sec that wc can even 
calculate a given connected correlation function g starting from the knowledge of the 
corresponding connected correlation function gj of the related Ising model. Hence, for 
example, in infinite dimension, for a correlation function of order two at zero external 
field, we have 

^ ^ I ^-P ({.TA\ (In: n: \ - In: \ In. \\ = fP. 



g^'\ii,i2) = J dV {{Jb}) ((a,;,a,;,) - (a,,) (a,,)) = ^^//(^i' ^2), (38) 
where g^fjjr{ii^ 12) is builded as in Eq. H24|) trough the correlation function g^p{i\^ 12) 

Qi (ii, 12) = (o-iiCTja)/ - {(^ii) I {(^12) I ■ (39) 

More in general, a connected correlation function of order k of the spin glass model is 
given by 

ff^'H^i,...,*^)^ j d7'({J,})(a,,...a,j(^) = .9f;_^(*i,...,zfe), (40) 

where (cTij ...aij.)^'^^ represents a connected correlation function of order k (see e.g. 
[TH]1. and g^^jf{ii, . . . ,ik) is given through 5^*^^ (ii, . . . , ik) with the analogous of Eqs. 
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I|23I24|I : g'f\ii, ■ ■ ■ ,ik) being the connected correlation function of the corresponding 
related Ising model 

gPi^,,...,^,)^{a,,...a,J^fK (41) 

In this work wc will not discuss the quadratic correlation functions 

dn{Jt}){{a.,...a0^^^r. (42) 

5. Disordered systems vs n interacting Ising like systems 

As is known, a way to carry out Eq. jSJ even with a fixed value of N (for shortness 
we will always omit the dependence on N), consists in applying the replica method 

/ dV HM) log {Z ({J,})) = hm = lim (43) 

where we have introduced the notation 

/ dV {{Jb}) Z {{J,}r . (44) 



= cxp -(Jb-Jo)V>/' , (45) 



Unlike Z", Z'") represents the partition function of n interacting Ising like systems, 
i.e., a system with n families of Ising spins cr'^', . . . , cr'"' interacting each one through 
the usual two-body Ising like interaction and, among themselves, through a non 
quadratic additional interaction. For example for a Gaussian measure 

d^i ( Jb) _ l_ 

it is easy to recognize that 

i^b ^■■■,Cb t 

where the replica indices a and /3 run over 1, . . . ,n. A part from constant factors, 
the case n = 2 of Eq. H46|l corresponds to the Ashkin- Teller model, for which in two 
dimensions some exact results are known |17 | -|19 ) . In particular, it is immediate to 
solve this model when one chooses Jq = 0, corresponding to the symmetric Gaussian 
measure. In fact, in this simpler case in Eq. H46() we are left with the only Ising 
variable Sf, = o'^^^aj^^^ and we are so reduced to a pure Ising sum as follows 

E /(•)= E /(■)' (47) 

E /({-I''-r})-2ir'E/(N})' (48) 

|r| being the number of bonds. Note however that for n > 2 the model Z^") is no 
more soluble even for a symmetric measure. This happens due to the impossibility 
to redefine suitable Ising variables, an aspect related to the concept of frustration. 
On the other hand, even in the general case, it is always possible to find Ising-like 
contributions in Z'"'. As we shall see later, if the dimension D is sufficiently high, 
above the critical temperature the Ising-like contributions become the leading part of 
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It is interesting to note here another aspect of the case n = 2. It has been well 
established that the Ashkin- Teller model presents a non universal behavior when the 
sign in front of the quartic coupling in Eq. H46() is reversed. In that case in fact, 
the specific heat diverges with a power law whose exponent is a continuous function of 
the quartic coupling. On the other hand such a non universal behavior cannot appear 
in the "disordered" system Z'^^\ simply because the measure (|45|l is definite only for 
real values of the disorder parameter J. This observation enforce then the idea that 
at least some features of a spin glass could be effectively described through a suitable 
Ising model. 

6. Representation of Z({Jh}) as sum over closed paths 

Let us consider a generic Ising model at zero external field with given couplings { Jb} 
defined over some set of links F. Note that, since the couplings arc arbitrary, what 
we will say will be valid, in particular, for the related Ising model. It is convenient to 
introduce the symbol 

Kb = pJb. (49) 
For the partition function it holds the so called "high temperature" expansion 

Z {{Jb}) = n cosh (Kb) 51 n (1 + t^i^li (^^^b)) ■ (50) 

As is known the terms obtained by expansion of the product Hber (1 + '^^fc tanh (Kb)), 
with k bonds proportional to ab^ db^ ■ ■ ■ '^bk j contribute to the sum over the spins only 
if the set 7 = {61, 62, ■ • ■ , ^fc} constitutes a closed multipoligon over F for open or 
periodic boundary conditions, and a collection of multipoligons and paths, whose end- 
points belong to the boundary of F, for closed conditions (when all the spins on the 
boundary arc fixed to be +1 or -1) (e.g. see [201); in such cases ab^^a-b^ ■ ■ -db^ = 1 so 
that Eq. H50|l becomes 

Z {{Jb}) = 2^ n cosh {Kb) n t^i^^ (-^f) ' (51) 

where the sum runs over all the above mentioned paths 7. Note that in the case 
tanh (iff,) = 0, the sum over the paths gives 1, {i.e. the contribution with zero paths 
must be included). Given the set of couplings { Jb}, in the thermodynamic limit, the 
series in the r.h.s. of Eq. (|51|) . normalized to iV, will be convergent for values of 
the parameters Zb = tanh(/3Jb) sufhciently small, i.e., inside a suitable set S! whose 
boundary corresponds to a critical surface S of the quenched Ising model. Since in 
Eq. (|51|l we have a power series, in Zb, with positive coefficients, it turns out that 
as anticipated, is a convex set. 

Let us consider the important case of a homogeneous model: Kb = K. In this 
case, Eq. ()51|l can be written in the form of a power series in the single parameter 
tanh(ii:) 

N 

Z = 2^ coshl^l {K) Ci tanh' (A') , (52) 

where Ci is the number of paths of length /. Note that, for fixed /, C; is a function 
of the size N and Ci oc N, so that in the thermodynamic limit the useful quantity 
is Q = limjv^oo Ci/N. Whereas Ci represents the total number of paths of length / 
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in a given lattice of finite size iV, c; represents the total number of paths of length 
I passing trough a single site in an infinite system; c; is a growing function of I and 
for large I this growth is known to be exponential, both for lattice systems and for 
tree-like structures, see respectively Refs. [21] and EH and references therein. As 
Ci is positive, for positive values of A', the series will be convergent and absolutely 
convergent if tanh(if) < w^^^ where 
1 



w 



- lim c/ , (53) 



F 

whereas for negative values of if, according to the Leibniz criterion, the series will be 
convergent if w^j^p < tanh{K), where is defined as the largest (in modulo) value 
of tanh(_ftr) such that 

lim Qtanh'(/v) = 0. (54) 

/ — >oo 

Hence, for J > 0, in the thermodynamic limit the series ci tanh^(/?J) exists for any 

( F) ( F) 

/3 < /3r % where the inverse of the critical temperature, is solution of the equation 
Wp^ = tanh(/3c^'' J), and analogously, for J < 0, the series will be convergent for any 
P < f3c^^\ where j3c^^\ the antiferromagnetic critical temperature, is solution of the 
equation w^^\, = tanh(/3c"^^^ J) . In the thermodynamic limit, given a coupling J, for 
any value of (3 lesser than the critical one, it is possible to define the typical length of 
the path Fby looking at the distribution pi = ci tanh(/3J)'/(^;, Q' tanh(/3J)' ). Note 
that, as — > P~ , the distribution pi becomes infinitely flat. Similarly, in the most 
general case, in which at any bond b we can have a different coupling Jb, the partition 
function can be rewritten as 

N 

Z{{Jb})^2^l[cosh{Kt) ^ C{{k})l[{tanhiKb)t , (55) 
ber {h=o} b&i 

where C{{lh\) is the number of paths having {/t,} units along the bonds {&}, with 
lb = 0, 1. Even in this case, from the knowledge of C{{lb]), it is possible to calculate 
the typical length per site li, along the bond b. As in the homogeneous case, also in the 
general case the critical temperature of the system is determined by the asymptotic 
behavior of the rcscalcd coefficients c({/b}) = limAr^oo C{{lb])/N . By rc-phrasing in 
terms of paths, in general we have that: 

The critical behavior of the system is determined by the paths of arbitrarily large 

length 

The above observation will be crucial when F is a generic graph infinite dimensional 
only in the large sense. 

7. Averaging over the disorder 

Let us now introduce the universal function P which takes into account the non trivial 
part of the high temperature expansion 

F({z,})^^n^^' (56) 

and let us average P over the quenched couplings (the disorder) 

({F,^^)}) ^ jdV ({ J4) F({tanh(if,)}) , (57) 
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a 

Figure 2. Example of one single path 7 on a 2-dimensional lattice. In this ease 
7 is constituted by three sconnected parts: a, b and c. Note that there is no 
overlapping between bonds. In the part c of the path we have drawn arrows to 
stress that this part of the path is constituted by only one single connected part, 
the path c being connected. 



where we have introduced 

= j d^ibtemhiKb). (58) 

From the product nature of the distribution dV {{Jt}), Eq. lO, it is immediate to see 
that P'^) is given in terms of the function P through 

P(i)(l^(i)})=p(l^a)|)^^j-j^(i). (59) 

7 be7 

Later, to evaluate the free energy we wih need to consider also the averages of 
Pf" for n e N 

[{fI'\ ^JdV ({J,}) PF ({tanh(7^,)}) , (60) 

where for m = 1, . . . , n we have introduced 

F^""^ = J dfib (tanh(A-,))'" . (61) 

We note that, according to Eqs. H13I14|I . unlike P({tanh(A'f,)}), the function 
-f ({^& is the non trivial part of the high temperature expansion of the related 

Ising model with couplings {Pj^"^^}. 

Let us now generalize Eq. H59(l to P'^"^ From Eqs. (|56() and H60() we see that 
for n integer we can calculate by summing over n replicas of paths 71, . . . ,7„, 
specifying for any of their bonds how many overlaps are there with all the other paths 
(see Fig. EJ. We arrive then at the following expression 

P^")= J2 jdVilh}) n tanh"(i^,)x 

n n tanh^"-!) (/n) X 

(n) benjLj_,_,;^7i\7ii 



n n tanh^"-^) (A', 

(n,J2) foen^^i ,^7A(7nU7i2) 
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Figure 3. A contribution to the summation of Eq. 1621 with n = 4. Here we 
have: a bond with overlap of order 3, 62 = 71 H 73 D 74; two bonds with overlap 
of order 2, bi U 63 = 71 D 73; and all the other bonds with no overlap (order 1). 
Note that the paths 71 and 72 intersect each other as geometrical objects but not 
as sets (for definition, a path 7 is the union of its bonds). The same observation 
holds for the paths 73 \ (f)i U 62 U 63) and 74 \ f'2- 



W W tanh(iffc), (62) 

(i^) he7i„\(U!^>„7!) 

where, in the product 11(11 12 ik)' ^^'^ indeces 11,12, ■■■ ,ik run over the n\/ ((n — fc)!A:!) 
combmations to arrange k numbers from the integers 1, . . . ,n, and with the symbol 
(ii,i2) we mean the couple «i,i2 with ii ^ 12 and similarly for (ii,i2, ■ ■ - ik)- From 
Eq. by using Eq. ((T)) and the definitions (jSTJ, we arrive at 

p'-'- E n ^"'x 

7i,...,7„ 6en;'^i7i 

n n p^'' X 

(ji) ben^^i ,^.^7,\7ii 

n n F("-^x...x 

(n,«2) ben;Li,i^ii,,27A(7iiU7i2) 

n n pi^'- (63) 

(ir.) h67i„\(Ui^.„7!) 

8. Free energy 
From Eq. (jSlfl we have 

dV {{Jb])log{Z {{Jb})) = J dV{{Jb})log |^2^ncosh(iffc)^ 

+ fdV {{Jb}) log I 5] n tanh(/u) , (64) 
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from which, by using Eqs. (0) and ((JJ) in the first term of the r.h.s., we get 

-I3F^N log(2) + ^ y" d/ifc log (cosh(A-b)) + 0, (65) 



where the non trivial part cj) is given by 

<j,= I dV {{.h}) log n tanh(/ifc) . (66) 

With the symbol 0/ ({4^^}) 

we will mean the non trivial part of the free energy of 

the related Ising model 

0,({4^)})^log(p,({4^'})). (67) 

The densities of (j) and 4>i will be indicated as Lp e Lpi, respectively. The free energy 
term will be obtained in terms of P("\ Eq. I)60|l . via the replica method with the 
analog of Eq. 

p(") _ 1 

6 = lim . (68) 

n-»o n 

9. Centered measure 

In this section we will consider the case of a centered measure: 

diibiJb)Jb = 0, (69) 
which in particular includes the symmetric case: 

dfJ-bi^Jb) = dfib{Jb)- (70) 
For a centered measure we have 

p{2m+l) ^ Q ^^^^ 

As a consequence, in the high temperature expansion, we find the following structure 
= 1, (72) 

p''' = E n -En Pb'' = pi{Fi''}), (73) 

71=72 bG7i 7 be-f 

E n 

7i.--->72,i6<f2n 66nfZ!i7i 

n n Ff"-^x...x 

(»i,»2) ften^^i .^7A(7,iU7,2) 

n n (^4) 

(z„„i,i„) &G7i„_in7i„\(u,^,„_j,,„7,) 

where (o2n is the set that constrains the sum in the r.h.s. of Eq. H74() to be restricted 
to combinations of paths for which any bond b may have only an even number of 
overlaps with any other bond. 
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Unlike the case of Eq. H7,S|I . corresponding to the centered Ashkin- Teller case, we 
cannot express the general 2n-th term of Eq. (|74|l in terms of an Ising like suitable 
sum as 

PUpI""^})- Nevertheless, for any n, m the r.h.s. of Eq. ^ we recognize 
Ising like contributions such as 

7 be-y 71,72:71072=0 be7i 6672 

7 bey 71,72:71072=0^671 bey2 

7l,72,73:7in7j=0,i7^j !)67i &S72 ^€73 

and in general we have 

0<7ni <m2 . . .<m„ :mi + . . .m„ — n 

where the coefficient C"^^"^ (2mi, . . . , 2mn) is given by 

(2mi, . . . , 2m„) = ^j^"'" , ^ ^ , (78) 

where n' is the number of different values of the m's in the sequence mi, . . . , m„, and 
the numbers gi, . . . , g^/ take into account of the degeneracy of the values mi, . . . , m„: 

„ _ / X]/=l '^mp,mi ! mp 7^ 

1 1, mp = 0, ^^^^ 



with p = 1, . . . , n'. The number 3 which appears in front of the second term of Eq. (|75|l 
counts the number of possible ways to pair (two to two) 4 paths. Similarly in r.h.s. of 
Eq. H76|l the numbers 15 in front of the second and third terms take into account of the 
possible ways to pair 4 paths and 6 paths, respectively, from a set of 6 paths. Finally, 
in Eq. H78|l we have written the general form of these combinatorial coefficients. 
However, as we will see, these combinatorial coefficients are of no importance in the 
thermodynamic limit. 

In Eqs. (|75I77(I we have explicitly written all the contributions in which any 
path 7 coincides with other paths an even number of times and we have then re- 
named the remaining different paths as similarly done in the third member of Eq. 
H73|l : 7 = 7i =72. AH these terms, up to the constrain that the re-defined paths 
7's cannot have common segments, are Ising like terms, whereas in Eqs. H75I77|I 
we have left dots, . . ., for indicating the other contributions which come from more 
complicated combinations of paths in which two or more paths 7 overlap by themself 
only partially {i.e. not entirely two to two) as shown for example in Fig. 0] These 
latter contributions cannot be expressed as Ising like terms. 



9.1. Approximation in high dimensions 

Up to now we have not yet introduced any approximation in our scheme. 
Unfortunately, the equation (|77|) for p(2") presents elements of intractability which 
avoid an exact calculation even in two dimensions where, on the contrary, the exact 
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Figure 4. A chain of three connected planar paths encapsulated in a larger 
planar path. The four paths overlap each other only partially. The paths in the 
figure arc slightly shifted for visual convenience. 



solution is possible for the pure Ising model. These difficulties are of two kind. First, 
in the Ising like terms, the re-named n paths 7 cannot have common segments or, in 
other words, we have to sum over the ensemble of n non overlapping random walks. 
Second, we do not have any knowledge about the non Ising like terms. 

Nevertheless, wc can show that, as the dimension D of the system grows, 
neglecting either the constrain for the non overlapping random walks and the non Ising 
like terms, implies a smaller and smaller error in evaluating the general expression (|77|l 
which, in the limit D — > cx), becomes exact. 

In the following we will make use of the fact that for finite and positive values 
of m we have 0(F(|^"^) - ©(F^^^)™. In particular, for the ±J distribution we have 
p{2m) _ ^p(2)-^rn^ anticipate that, in the limit of infinite dimension, any finite 
difference between F^^^™^ and (i^^^'')'" becomes irrelevant for the mapping to be valid. 
For r, to be specific, we will consider a D-dimensional hypercube lattice. However, 
as will be evident later, similar arguments can be repeated for any systems of links F 
infinite dimensional in the large sense (see Sec. 3), as happens, but not only, in the 
case of generalized tree-like structures. 

For the moment, let us consider for simplicity only planar connected paths. Let 
us start with n — A. We have to sum over all contributions that belong to the set 
<£i. From the set of all the possible planar connected paths, wc have to draw out four 
arbitrary paths and combine them in any possible way compatible with the set (S4. 
We find useful to decompose this set as follows 



where (04^4 is the set of 4 coinciding paths, i.e. the set of all one replica paths, see Fig. 
121 (02,4 is the set of all paths overlapping in couples (two to two), see Fig. |H| and ^4 
is the rest, i.e. the set of all paths in (04 which overlap each other only partially, see 



Wc now observe that in general, if D is the dimension, and _D > 2, at any point 
of our lattice system we have D orthogonal planes where anyone of the 4 planar paths 
can be arranged. On the other hand, whereas any given element of (04^4 or ^4 lies 
necessarily on one single plane, we can generate an element of (92,4 by arranging two 
arbitrary paths, even in two separated planes, for a total of D{D — l)+D combinations. 
Therefore, we have 



£•4 ^ (?4.4Ufr2,4U^4, 



(80) 



Fig. H 




(81) 
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Note also that, up to the constram 71 H 72 = 0, the set (02,4 corresponds to the range 
of the summation in the second term of Eq. H75|l . 

Analogously, for n = 6, we have to consider the following decomposition 

J'6 = <f6.6UA,6UJ'2,6U^6, (82) 

and we see again that the elements of the set (02,6 can be arranged in D{D — 1){D — 
2) + 0{D{D — 1)) ways, whereas the elements of the other sets can be arranged at 
most in D{D — 1) combinations. Therefore we have 



Note that, up to the constrain for the non overlapping of the re-named paths, the set 
^2,6 corresponds to the range of the summation in the third term of Eq. (|76|l . 

We can repeat the above argument for any positive integer n. We always arrive at 
the conclusion that the Ising like contributions obtained by pairing in all the possible 
ways two paths are the leading contributions of the set as follows 

K2„| - K2.2„| ^ / 1 \ (34^ 



Up to the constrain that the re-named paths must not overlap among themselves, 
the set (o2,2n corresponds to the range of the summation over paths in Eq. (|77|l . On 
the other hand, by applying the same argument of the dimensionality, we see that 
neglecting the constrain implies a vanishing error for large D, being 

K2,2n| — X]7i,...,7„:7,n7j=0,i#j 1 1 



Finally, it is not difficult to convince oneself that the a very similar argument 
can be applied also for general paths, non connected and non planar. In conclusion, 
taking into account Eqs. and (P^ . the general expression (|77|. and the fact that 
0(F(2™)) _ 0(F(2))™, we arrive at 

Furthermore, since for » iP^J^)-^ the typical length 1{P) goes to zero, we 
can improve the description of the above error including an unknown bounded function 
g{x) of order 1 which goes to zero for /3^^ >> (/Jc^"*)^^ and takes into account the 
fact that in such a limit Eq. H86|l becomes a trivial identity, being both and 

P" ^{-Ffl^''}^ equal to 1 (as we shall see shortly the coefficient in front of P" does not 

play any role in the thermodynamic limit). Therefore we can replace Eq. H86|l with 
the finer one 

^■■"-M^"(<^"))^o( ^""'t:jf'-'' )- m 

By exploiting this approximation, we are now able to calculate the free energy term 
(j) via the replica formula of Eq. H68|l or, more precisely, via the following limit 

p{2n) _ I 

(h = lim . 

n— 2n 

By using the relation 

{2n)\ = 2"n!(2n- 1)!!, 
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and the Gamma function property 

r (n + i) 

(2n~i)!!^ r(i) 

form Eq. (|57jl we arrive at 



P(-) = ^^jTf^2-p- + O ) . (89) 



Finally, by taking the derivative with respect to n we get 



p{2n) _ 2 I 



= lim : 1 = 1 log ( p ( \ ^ j ^ 



liog(2) + ird,!2£<2, (90) 

^ ^ Jo i^e*^ 



where we have used 



Equation ()90|1 can also be written as 

Vf })+0(^)+^iog(2) + ;: / dt:^. (92) 



Finally, we note that in the thermodynamic limit the last two terms in the above 
equation are of no importance, so that, for N large and by using Eq. H65() . for the free 
energy F we find the following expression 

-PF=J2[ dM6log(2cosh(i^,)) + ^cf>j ({Pf '} 

+ o(^^y D>2, (93) 

Pf ^ = J dnbtanh^{Kb), beT. (94) 

Equation ()93I94|) has been derived starting from the high temperature expansion of 
the related Ising model, therefore, taking into account that ^/ is convex, it holds for 
any temperature /3~^ such that / dfXb tanh^(/3J{,) < Wb- Equation H93I94|I tells us 
that for such temperatures, the free energy of an Ising spin glass model defined over a 
D-dimensional set of links F, with D large and in the presence of a centered disorder, 
can be effectively expressed in terms of the free energy of an Ising model defined over 
the same set F in which the parameters of the high temperature expansion are replaced 
by the following effective substitution 

Zb = tanh{Kb) Pf ^ = J di^ib tanh^ {Kb), beT. (95) 

In particular, if the set of equations J djib tanh^(/3Jb) = Wb admit a solution for some 
j3c, ^ill be the critical temperature of a P — SG transition. 
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Figure 5. Schematic example of four completely overlapping planar paths. The 
paths in the figure are slightly shifted for visual convenience. 



Figure 6. Schematic example of two couples of two completely overlapping 
planar paths. The paths in the figure are slightly shifted for visual convenience. 



10. Generalization to non centered measures 

If ^tjI^^+i) is more 0, in constructing P^^n)^ besides the terms in the set S'2n, in which 
appear only terms involving f^^^™-*^ rn = 1^ . . . ^n, we have to consider all the other 
terms belonging to the complement of S'2n, ^im *-e-, all the families of contributions 
in which among the 2n paths, there is at least one path in which one or more bonds 
do not overlap any other bond of the other paths or, if there is an overlap, it is an odd 
overlap, i.e., an odd number of paths overlap over the same bond. Even in this case 
we can single out the Ising like terms either with even or odd overlaps of the bonds. 
To this aim we find useful to decompose the whole set of 2n paths as: 

U #2„ = ^0 U . . . U J^fe U . . . U ^„ U ^„ (96) 

where, for each fc, is the set in which any element is constituted by 2n paths in 
which 2fc of them overlap each other only an even number of times and the remaining 
2?! — Ik paths overlap each other only an odd number of times, whereas Sin is the set 
of all the possible non Ising like terms. 

10.1. Approximation in high dimensions 

When the dimension D is sufficiently high, in each subset we can repeat the 
same argument as in the previous section as follows. Among the 2k paths (only even 
overlaps), the most important terms are those with minimum non zero overlap, i.e., 
those terms which give a factor oc P'^HFj^^^}), and, similarly, among the remaining 
2n — 2k paths (only odd overlaps), the most important terms are those with no overlap, 
i.e., those terms which give a factor oc p2"-2fe(|^a)|)^ ^^^^ ^^lat the proof for the 
second statement runs as in subsection H9.1|l . the only difference now is that it works for 
paths having odd overlaps, so that, in this case, the minimum number of overlapping 
paths is 1, as opposed to 2 in subsection (|9.1|) . and, as a consequence, the relative 
error in neglecting the other terms goes as 1/(1? ~ 2n+ 1), instead of 1/(Z) — n + 1). 
Furthermore, the minimal dimension to apply the argument for paths with only odd 
overlaps is not D > 2, but D > 1. 
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For what concerns the set we can neglect these non Ising like contributions 
also by applying almost the same argument of the subsection l|9.1() . The only difference 
now is that these non Ising like terms involve also paths having both even and odd 
overlaps. 

In conclusion, wc arrive at the following expression for p(^") 

fc=o ■ \ / \ / 

where the factor (2A:)!/(2''"fc!), as before, takes into account the number of ways to pair 
2k paths, whereas the binomial coefficient takes into account the number of ways to 
choose 2fc paths from a total of 2n. Yet, as in the previous section, these combinatorial 
coefficients are of no importance in the thermodynamic limit, where, both the factors, 
P({Pj|^'}) and P({fI^^}), grow exponentially in N, so that, in this limit, only the 

two leading terms, P"'{{fI'^^}) and P'^^{{Flp}), are important and, by taking the 
logarithm of this sum in Eq. H88|l , the final formulas given in Sec. 4 follow. 

11. General graphs 

In the previous sections we have proved the mapping by providing a simple argument 
which applies when the dimension D is sufficiently high. Roughly speaking, the key 
point is that, if D is the number of axis passing through a site, by choosing at random 
two of them (or replica), the probability that they coincide (or overlap), goes to zero 
as 1/D. For pedagogical reasons we have used this argument by referring to models 
whose set of links F is defined over a D dimensional hypercube lattice A, where D is 
related to the number of first neighbors, 2D. In these models it is easier to visualize 
the axis. However, as we will see soon, little changes are involved in the proof if we 
consider a set F infinite dimensional in the large sense. Even for these systems, we 
can always find the analogous of the number D of independent axis per site and look 
at the situation when D ^ oo. 

Let us see now, more specifically, an Ising model defined over a Cayley tree of 
coordination number q = k + 1. It is constructed as follows: one starts from a vertex 
root '0' and adds q points all connected to 0. This set of q points represents the first 
shell of q sites. The second shell is instead obtained connecting each one of these sites 
to new k = q — I points, and so on for the successive shells {k is then the branching 
factor) . 

This lattice is a tree and, therefore, if we do not close in some way the boundaries, 
or if we do not broke in some minimal way the symmetry up-down, as happens on 
a hypercube lattice, the non trivial part of the free energy ipj is zero, and no phase 
transition is possible. Note that, since a Cayley tree depends heavily and non trivially 
on the boundary conditions, what is usually studied is a Bethe lattice, consisting in 
a subset of the Cayley tree infinitely far from the boundary. Yet, even with such a 
definition, it is known that, unlike the Ising model over the Bethe lattice, the spin 
glass model over the Bethe lattice remains still dependent on the kind of the boundary 
coinditions imposed on the sites of the outer part of the Cayley tree. It is known in 
particular that, even though in this model there is a singularity in the free energy, 
such a singularity implies a "true spin glass phase transition" (i.e., with the same 
qualitative picture as the replica symmetry breaking in the SK model ^) only under 
certain boundary conditions However, such a distinction for us is not important; 
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our aim here is limited to find only the singularities of the model, regardless on the 
fact that such singularities may or may not imply a true spin glass phase transition. 

Let us consider a centered measure and start from the general exact representation 
H74|l . As we will see, unlike hypercube lattices, if we analyze the Ising- like contributions 
for a Bethe lattice, we do not reach exactly the Eqs. H75I77|I . but a slight modified 
version of them, which, however, in the thermodynamic limit, coincide. Let us look at 
the set of all the possible paths starting from the root and arriving at some infinitely 
far boundary. Note that, in a Bethe lattice, regular or not, such a set coincides with 
the set of all possible paths. Furhtermore, as will be clear soon, in a Bethe lattice, 
and more in general in a tree, due to the absence of loops, the effectes of non-Ising 
like terms are irrelevant, being limited to an overlapping of a finite number of bonds. 

It is immediate to see that the term P^^^^ remains formally as Eq. (|73|l (see Fig. 
Ej). Let us calculate P^'^\ Due to the absence of loops, if 4 replica paths starting from 
coincide for I bonds, after a splitting at the bond ? + 1 in 2 branchs of 2 coinciding 
paths, the branchs will not overlap each other anymore (sec Fig. |HJ). Therefore, taking 
into account that the total number of paths of length / starting from 0, is gfc'"^, the 
equivalent of Eq. (|75|) for a Bethe lattice with n shells becomes 

n I 

1=0 r=l 71.72 be7i\u;.^i6^ fce72\u;,^ib^ 

where {br}r=i, is a sequence of n successive bonds starting from 0. On the other 
hand, in the thermodynamic limit, n oo, the constrains \ ul^i br in Eq. (|98|l 
become negligible and we get 

p(^)=3gf:fc'n<v^({^^f}). (99) 

1=0 r=l 

Similar expressions can be derived for any n. Even if these expressions become quite 
involved for increasing values of n, in the thermodynamic limit we are always left with 
a general form of the type 

where Bn is a suitable constant which takes into account the combinatorics, Eq. l(7^ . 

and integers powers of the sums (zX^z^o ^' nt=i -^b™^ for 4 < m < n. Since both 
of these factors do not grow with the size of the system, as for the previous case of 
subsection 9.1, we arrive at the same conclusion of Eq. (|93|) (with £> = oo in this 
case). 

A similar argument can be repeated for tree-like structures and for graphs in 
which there is at most a finite number of loops per vertex. Unlike the pure Cayley 
tree now, we have to consider also a finite number of closed paths (see Fig. EJ, but 
the construction of the terms p(^") runs in the same way. Finally, we see that, as 
anticipated in Sec. 3, we can consider even more complex non tree-like structures in 
which the number of loops per vertex is not finite; the only condition we need being 
that the paths, two at two, share at most a finite number of bonds. The only problem 
with such complex structures is that the related Ising models are hardly solvable. An 
exception to this is provided by the fully connected graph, i.e., the SK model. 
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Figure 7. An example of a Bethe lattice with coordination number g = 3 
(fc = 2). The path of greater thickness represents two completely overlapping 
trajectories 7 = 7i =72 going toward the boundary of the lattice. 




Figure 8. Bethe lattice with tj = 3 (fc = 2) with four trajectories overlapping 
four times along the bonds with label "4" and overlapping 2 to 2 along the bonds 
with label "2". 

12. Applications to finite dimensional models 

We now apply the mapping to some finite dimensional cases and compare them with 
known results. Hereafter we will consider mainly homogeneous models, d/ib = d/i, 
so that, in particular, .FjJ^' = F^'^\ Furthermore in this section we will consider only 
centered measure and hypercube lattices. For the Ising like term, ipi{F^'^^), in one and 
two dimensions, we can use its analytical knowledge coming from the exact solution 
of the Ising model, whereas, in higher dimensions, we can use for it the knowledge 
coming from numerical methods. Clearly, since our approach takes into account only 
the leading term of a 1/D expansion, in low dimensional systems, we expect poor 
results, but it is however interesting to see even in these cases how our approach 
works. We will force the study even in the cases D = 1 and D = 2 which can be seen 
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Figure 9. A generalized tree-like structure with loops obtained starting from a 
a Bethe lattice having previously q = 3 (A: = 2). The drawn loops pass through 
the root vertex 0. The numbers over the loops are the corresponding lengths. 



as singular points where the mapping is not definite. 

12.1. One dimensional case 

For D = one neglects irrelevant boundary effects, there are no closed paths so 
that from a direct calculation one sees that Lp is exactly 0. The one dimensional case 
is therefore trivial. Nevertheless, it is worth to observe that, for the same reason, even 
ipi is exactly zero, so that for Z? = 1, the mapping, even if not definite, turns out to 
be exact. 

12.2. Two dimensional case 

For 07 in two dimensions we can use the Onsager's solution. If we indicate with 
Zfi = tanh(fc;i) and z„ = tanh(fc„) the horizontal and the vertical parameters of the 
model, Onsager's formula in the thermodynamic limit gives |25| 



In the case of an isotropic spin glass, we have d^h = dfiy, so that it is enough to 
consider the related Ising model with z = Zh = z^,, and Eq. (|101() simplifies in 





■2 



■1' 



2zhZy {zh COS {9.„) + z,u cos {Oh)) + -2zh cos {0h) - 2z^ cos {9y)]. (101) 




(102) 



where 



= Az 




(103) 
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As is known, (fi (z) is non analytic at k = ±1, for which the specific heat, as a 
consequence, has a logarithmic divergence in the temperature. In the high temperature 
expansion representation, the condition k ~ ±1 amounts to w^^' = ±('\/2 — 1). 
This implies that in D ~ 2 the pure Ising model presents a ferromagnetic or 
antiferromagnetic second order phase transition at the critical value kc^^ = 0.4407 
{Tc — 2.269|J|, if the Boltzmann's constant is taken as 1). From Eq. (|93|l we see 
therefore that our approach for a two dimensional Ising spin glass would predict a 
second order phase transition at (3c solution of 

j rf/xtanh2(/3(^G)j) = ±(V2- 1). (104) 

Of course, this equation may have solutions only if the r.h.s. is positive so that 
only a correspondence with a "ferromagnetic" -like transition is possible. In particular 
for the plus-minus distribution © centered in J and —J, a second order phase 
transition would take place if tanh^(/3c'^'^' J) = \/2 — 1, which has solution for 
ki^'^^ = (i'f'^^J = 0.7642 (ri^^' = 1.308J). This result is in contrast with the (by 
now) known fact from numerical simulations that, in two dimensions, the Ising spin 
glass has no phase transition at finite temperature )26| . As is evident from subsection 
9.1, D = 2 represents a singular case where the mapping is not definite since inside 
a bidimcnsional space there is only one plane. This explains why in our approach in 
two dimensions we find a finite critical temperature; for the two dimensional case, an 
effective mapping with a suitable Ising model is impossible. The fact that there is not 
a finite temperature phase transition, implies that in two dimensions the non Ising 
contributions not only are not negligible, but in the thermodynamic limit constitute 
the leading part of the high temperature expansion and hide the Ising like effects. 
Nevertheless, it is interesting to see that some features of this model may be explored 
even in our Ising-like approach. For example, as appears by comparing the plot of 
Fig. EH foi' the internal energy with the exact numerical results obtained in |2Z| , 
though the internal energy is wrong, from its slope, away from the critical point, 
we get a certain evaluation of the specific heat. Note also that, even though it is 
clearly wrong the result for the critical temperature, according to the general rule 
H34|l . in our approximation, the effect of the disorder has however decreased the critical 
temperature [Tc^'^^ = 1.308J) with respect to the critical temperature of the pure 
Ising model (T^ = 2.269 J). 

12.3. Three dimensional case and higher dimensions 

For D > 2 there is no analytical solution either for the pure Ising model and for 
the Ising spin glass model. Nevertheless, several very accurate numerical (and partly 
analytical) data are nowaday available for D = 3. From high and low temperature 
expansions and other numerical techniques it is known that in £) = 3, the pure 
Ising model has a critical point at ki'^ = 0.221 {Tc = 4.51154J) [SS], so that 
w'^^ = tanh(fcc^'') = 0.2180. Therefore, from the general rule (|19(l . we obtain that 
in D = 3 an Ising spin glass has a phase transition at a critical temperature /3c^''' 
such that 

dut&nh^iPi'^^U) = 0.2180. (105) 

In particular, for the plus-minus distribution, the above equation gives kc = 
pi^^h = 0.5062 (Ti'^^' = 1.975). On the other hand from gHj a critical temperature 
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Figure 10. Internal energy normalized per bond and per unit coupling, e = 
E/{JNi,), for the two dimensional plus-minus spin glass model. The calculation 
of e is reported in Appendix B. 



at fc^^^^ = 0.851 (ri"^^' = 1.175) is found. 

Similarly, in higher dimensions, we find more and more agreement between our 
scheme and known results |3D1 131L , the difference being in good accordance with 
our estimation of the relative error whose order is expected to be l/D. 

For sufficiently high, but finite dimensions, we can write a general formula. Let 
us consider an hypercubic D dimensional spin glass model with a disorder having 
variance 1/2D. According to Eqs. I|12I14|I . the related Ising model has the following 
Hamiltonian 



Hr 



ber 



J(7b, 



(106) 



so that, if D is large enough it has a critical temperature located a |22 

(107) 

we find that the critical 



Therefore, by expanding at the first order in l/D Eq. (|19|1 . 
temperature for the corresponding spin glass is given by 

1 



= 1 + 



(108) 



in accordance with the result of references |10) and |33| in which a perturbative 
expansion in powers oi l/D is performed. 



13. Applications to infinite dimensional models 

As D oo our scheme becomes exact. In particular, for the SK model, above the 
critical temperature it reproduces the mean field solution, even though the procedure 
is completely different from the replica and cavity methods Notice that, as we have 
just anticipated in the Sec. 4, our approach, despite of the fact that in Z? = oo becomes 
exact, gives access to the free energy only above the critical temperature and, an 
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analytic continuation below the critical temperature is not allowed. Note, in particular, 
that a forced use of the mapping to calculate the free energy at low temperature, would 
give a completely wrong result. In fact, given a system in D dimensions, it is easy 
to see that ipj{P)/(3 ^ for /3 ^ oo, so that, at zero temperature, from Eq. (|65|l 
applied to a plus-minus measure, it remains only Fq/N — Uq/N — —JD, whereas 
the Derrida lower bound for large D gives Uq/N ^ —J^2D log(2) However, as 
just pointed out in the subsection 4.3, we argue that an analytic continuation of other 
physical quantities, such as the crossover surfaces or the magnetizations, provide a 
certain effective approximation. 



13.1. Sherrington Kirkpatrick model 

The SK model corresponds to the spin glass over the fully connected graph and 
anyone of the N spins interacts with any other spin trough a random coupling Jf, 
whose probability distribution has homogeneous rescaled mean value and variance 
given respectively by 



dfi Jb = Ja/N, (109) 

dfi {Jb - Jo/Nf P/N. (110) 
From Eqs. (|13|) we see that the Hamiltonian of the related Ising model is 

Hi = -Y, J^'^ab = -J2 (111) 

where in the last expression we have rewritten Hj in the usual form as a sum over the 
couples of sites (i,j). As is well known, for this model, depending on the sign of the 
coupling J^^\ a ferromagnetic-paramagnetic or an antiferromagnetic-paramagnetic 
phase transition takes place at the same critical temperature given by 

/3«|jW|7V = l, (112) 

which for N large, in terms of the universal quantities Wpjj^p = ± tanh(/3c^'' | J*^-^) |) 
gives 

^FlF = ±i + 0(^) (113) 



'F/AF - ^ 

On the other hand by using Eqs. H109I110|) for TV large we have 



PjY ( 1 

d/i tanh^ (/3J,) = + O f _ ) , (114) 



and 



/ 



d/z tanh (/3,A) - ^ + O ■ (115) 

Therefore, from Eqs. (|19(l and (|2()|l . in the limit N ^ oo, we get the following spin 
glass (SG) and, depending on the sign of Jq, ferromagnetic (F) or antiferromagnctic 
{AF) phase boundaries 

/3f<=) J = 1 (116) 
/3fMf^)jo= ±1. (117) 
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Finally, according to Eq. by taking the envelope of the curves and H117|l 

we get the upper phase boundaries shown in Fig. 1111 In the same figure we report 
also the coexistence curves SG — F and SG — AF derived from the systems (jSH 
whose solution is given by 



f3J) =P\Jo\- (118) 

Note however that, unlike the upper phase boundaries, the coexistence curves arc not 
exact; they are only representative of the true coexistence curves. 

In infinite dimensions, our approach becomes exact only above the upper phase 
boundaries, where the high temperature part of the free energy, ip, becomes trivially 
0. On the other hand, an analytic continuation below these boundaries is not allowed. 
This fact can be understood considering that the thermodynamic of the related Ising 
model (|111() turns out to be well defined only above the critical temperature, which 
in the thermodynamic limit is infinite (/3c^^ — *■ 0); for a well defined thermodynamic 
the J*-^) in the Hamiltonian (|111|1 should be replaced by J'-^^ /N. Nevertheless, via 
analytic continuation, besides the extrapolation for the coexistence curves (|118|l . a 
simple estimation of the Edward Anderson parameters and of the magnetizations is 
possible. In Appendix C we show that in the SK model, the square root of Edward 
Anderson parameter and the magnetization, respectively indicated as msG and ttt-f, 
naturally emerge as effective fields in correspondence with the fields of the related 
Ising model. For Jq > 0, they are solution of the following mean field equations, 
respectively 



msG = tanh (^(/3J)'msGj , (119) 

and 

mp ^ tanh (/SJomp) . (120) 

Note however, that the above expressions are not exact. In particular, the slope of 
msG near the critical point is wrong. 

Similarly, when Jq < the square root of the Edward Anderson parameter and 
the magnetization, which now are described by two effective fields, to^"-' and m*-''-' 
related to the two sublattices a and b in which the initial lattice can be decomposed, 
are solutions of the following mean field systems, respectively 



m^g^ ~ — tanh {{PJ^'TfigQ 
m-gQ = — tanh ( {(3JYm"gQ ^ 

and 

fn'AF ^ ^ tanh {pjQm'A^r 



m 



AF 

(b) 
AF 



= — tanh 



(/3Jo 




(121) 



(122) 



13.2. SK generalized; random antiferromagnets 

One interesting generalization of the SK model, is a model defined over a lattice which 
can be decomposed in many, say p, sublattices eachone constituted of N sites and 
with the constrain that a spin over one sublatticc interacts only with the spins over 
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Jo/J 



Figure 11. Phase diagram of the Sherrington Kirkpatrick model. The upper 
Une, i. e. the envelope of the P — SG, the P — F and the P — AF lines, obtained 
from Eqs. 11161 and 11171 . is exact. The crossover curves SG — F and SG — AF 
are obtained from Eq. 11181 . 



the other sublattices. Such systems are particularly interesting as models of random 
antiferromagnets. The Hamiltonian is given by |3()| 

(123) 



H 



where the latin indeces i,j = 1,...,7V, label the sites in each sublattice, whereas 
the greek indeces u,i',= 1, . . . ,p, label the sublattices. A couphng J/^''^^ is then in 
correspondence with the bond b ~ (i, ^; j, ly) connecting the site i of the sublattice p 
with the site j of the sublattice v. Let us consider for simplicity a distribution with a 
homogeneous rescaled mean value and variance given respectively by 

Jo 



dp J^''^ 



dp 



Jn 



N{p-1)' 

N{p-iy 



(124) 
(125) 



N{p-l) 

The related Ising model has Hamiltonian 

Following |37] , the most general solution of the related Ising model must be found by 



introducing p effective fields my 



(1) 



,(p) 



satisfying the system of equations 



^tanh (-/3mf' 



N 



l^h 



(127) 



Note that the sign in front of Hj is, for convenience, reversed with respect to [HZI; we 
recall that the coupling j'^^ of the related Ising model is homogeneous, but its value 
can be arbitrary. Note also that, with respect to the convention adopted in the 
sign for the parameter Jq has been reversed. 

Linearizing Eq. (|127(l for small fields we arrive at the homogenous system 



where x is given by 

X = 



= —X TO 



(fc) 



; = i,. 



(128) 



(129) 
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The homogenous system H128|l has a non zero solution for values of x for which the 
matrix of the coefficients A has determinant zero. It is easy to see that this determinant 
is simple given by 

det{A) ^ {I ~ xf-\l + {p - l)x) , (130) 

therefore, the system will have a non zero solution for the critical values a: = 1 or 
X = —l/(p— 1), which using Eq. (|129|1 means 

-/3(^)j(^)iV= I (131) 

Recalling the explicit sign in front of the Hamiltonian l|126|l . we note that the first 
solution corresponds to an ordinary Ising ferromagnet with J^^^ > 0, whereas the 
second one corresponds to a generalized antiferromagnetic Ising model with J*-^^ < 
and with p sublattices. For N large, in terms of the universal quantities tanh(/3c'^'' J^-^^), 
Eq. (fnT|l gives 

^'aI = (133) 

which, according to the general rule (|19I20() and by using Eqs. (|124I125|I bring to 
the following spin glass {SG) and, depending on the sign of Jq, ferromagnetic (F) or 
antiferromagnetic (AF), phase boimdaries 

=1, (134) 

PPJo =1, Jo>0, (135) 

Pi^^^J„= -ip-l), Jo<0 (136) 

Finally, according to Eq. (|21|l . by taking the envelope of the curves (|134I136|I we get 
the upper phase boundaries shown in Fig. 1121 In the same figure we report also the 
coexistence curves SG — F and SG — AF derived from Eqs. (|35I36() whose equations 
are given by 

(/3j)'=/3Jo, Jo>0, (137) 
for the SG — F crossover, and 

(/3j)' = ^, Jo<0, (138) 

for the SG — AF crossover. 

Note that for p = 2 one recovers the standard SK model. 

13.3. Spin glass over Bethe lattice 

Let us consider a spin glass model defined over a Bethe lattice of coordination number 
q {i.e., branching number k ~ q — \, see definition in Sec. Ill|l . In the case of a 
homogeneous model, d^i, = d^, the related Ising model is the homogeneous Ising 
model over a regular Bethe lattice with coordination number k for which the exact 
solution is known This solution predicts a second order phase transition at a 

value of /3c^^ given by 
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Figure 12. Phase diagram of the generaUzed Sherrington Kirkpatrick model 
tT23l with p = 5. The dia gram can be compared with that reported in )36l . 
The upper Une, i.e. the envelope of the P — SG, the P — F and the P — AF 
lines, obtained from Eqs. 113411361 . is exact. The crossover curves SG — F and 
SG - AF are obtained from Eqs. IT?1 and 



In terms of the universal quantity Wp ~ tanh(/3c -Z^^^), Eq. (|139|l reads 

w^P = tanh J(^)) = (140) 

where we have made use of the hyperbohc identity 

tanhQlog(r)^ (141) 

Therefore, by using the general rule H19I20|I wc find that a spin glass {SG) and a 
disordered ferromagnetic (F) transitions take place at (3^^'^^ and Pc^\ respectively 
solutions of the two following equations 

Jd^l tanh^ = (142) 

j dii i'An\i{pf^Jf^ ^^'^^^ 

Finally, by using Eq. H21() and Eqs. (|35I36() . the upper phase boundary and the 
crossover surfaces follow. 



14. Conclusions 



In the framework of the high temperature expansion, we have derived a general 
mapping between an Ising spin glass model and a related Ising one, Eqs. H19I24|I and 
H25I30|I . The mapping is definite above the critical temperature and becomes exat in 
all the paramagnetic region when the dimension 13 — > cx) in the strict sense, whereas, 
more in general, becomes exact in the limit /? — > when the dimension D ^ oo 
in the large sense, where now D, roughly speaking, is the number of independent 
paths per site (see Sec. 3). The mapping can be applied to find exactly the upper 
phase boundaries of a spin glass model if the critical temperature of the related Ising 
model is known. Furthermore, by analytic continutation, the mapping provides even 
information such as the crossover surfaces, the magnetizations, and the correlation 
functions. Even if this further application is not exact, we argue that its use can turn 
out to be quite effective as a first insight to the physics of very complex models. 



Ising spin glass models versus Ising models I 



33 



We have applied the mapping to several known models at finite and infinite 
dimensions. In particular, the last application of a spin glass model defined over a 
Bethe lattice suggests a comment. We observe, in fact, that the relations for the Bethe 
lattice, Eqs. H142I143(I . have been derived by many authors in several contexts and by 
using different procedures, see for example EH -03- However, it is quite impressive 
the easiness by which we have derived these formulas (and similarly those for the SK 
models); simply starting from the solution of the related Ising model. For example, 
within the limits that our approach concerns only the singularities of the free energy, 
we had not to be worried about the rigorous definition of the Bethe lattice. As is 
known, in fact, the so called Baxter exact solution of the Bethe lattice concerns, more 
precisely, the inner part of an infinite Cayley tree far, in a sense, from the boundary 
in which a finite fraction of the total number of spins resides (see Sec. If one does 
not exclude in some way these boundaries, the model turns out to be greatly sensitive 
to the boundary conditions and difficult to treat. Another way to avoid this drawback 
consists in considering for example an ensemble of uncorrclatcd random graphs in 
which any site has a mean connectivity equal to q. We stress however that, as far as 
one considers only the singularities of the free energy, whatever the lattice and the 
system of links F over which we build a spin glass model could be complicated, in our 
approach, all the mathematical and technical difficulties, even those connected to the 
delicate questions concerning boundary conditions, if any, are completely reduced at 
the level of the corresponding related Ising model, which, as a non random model, 
is remarkable very simpler to solve. It is also clear that this difference becomes 
particularly important if the given model has many degrees of freedom (mechanical or 
disorder parameters). In our approach, once the phase boundary of an Ising model in 
infinite dimension is known, this solution can be immediately applied to find the upper 
critical boundary of the corresponding spin glass model for which the Ising one turns 
out to be the related Ising model according to the definition given in Sec. 4. Hence, 
for example, since the Ising model over a Bethe lattice dos not depend on the chosen 
boundary conditions, the upper critical surface of an Ising spin glass model over the 
Bethe lattice will be the same for any chosen boundary condition as well, regardeless 
of the fact that the non paramagnetic regions can be different (see comment in Sec. 
11)- 

We want to poin out also the numerical advantage that our mapping implies. 
In fact, even if the analytical knowledge of the critical temperature of some Ising 
model in not known, an its numerical evaluation turns out to be hugely easier than 
a direct numerical evaluation for the corresponding Ising spin glass model. Once a 
numerical estimation for the critical temperature l//3c^^ is known, the mapping returns 
immediatly the upper critical surface of the spin glass. 

The generalization of the mapping to include also a randomness of the set of links, 
necessary to study random models on random graphs (see, e.g., @31)i will be presented 
in a forthcoming work (part II) . Finally, we observe that the simple argument we have 
used to prove our mapping, "random Ising system " — > "non random Ising system" , 
which is exact in infinite dimensions, does not seem to be peculiar of the Ising model, so 
that, even for more general models, a mapping at high temperature "random system" 
— !■ "non random system" appears possible. 
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Appendix A. Free energy in the presence of external fields 

The result of Sec. 9 can be formally generalized for arbitrarly external fields {hi}. In 
this case, Eq. (|5n|l is to be modified as follows 

N 

Z{{Jt};{h,}) = n cosh (/u) J] cosh (if,) 

6er 4=1 

N 

X ^ [| (1 + ffb tanh (Kb)) n (1 + tanh (H,)) , (A.l) 
{rTi}ber i=i 

where we have introduced the symbol 

H, = ph,. (A.2) 

Correspondingly, the high temperature expansion of Z will be a series over arbitrary 
combinations of closed and open paths. Given one of these generalized paths, the 
weight of a bond 6, as before, is provided by tanh(irh), whereas the weight of two 
extremal points i and j of an open path are provided by tanh(ifi) and tanh(_ffj), 
respectively. Taking into account that the factors tanh(/fi)'s are not affected by 
the integration over the Jf,'s, the generalization of Eq. H63|l to generic external field 

iv) 

follows straightforward, where, along with the bond-weig htsF^''^^ now appear also the 
point-weights tanh'^(_ffi), where p and q count, respectively, the number of overlapping 
bonds and points, among the n replica generalized paths. On the other hand, it is not 
difficult to see that we can apply the argument on the dimensionality even to open 
paths and to combinations of these with closed paths so that even in the presence of 
an external field the mapping of Eqs. (|22I24() follows. The problem here is that this 
argument is only formal. In fact, if a nonzero magnetic field is present, with respect 
to the parameters Zb = tanh(_ft'f,), the high temperature expansion in general can have 
a zero radius of convergence when D = oo, and this is the case for the related Ising 
model; for hi ^ Q one has in general a non zero mean magnetization, so that the 
density energy of the related Ising model goes to infinite as _D ^ oo. 

Appendix B. Internal energy for the D = 2 spin glass 

In this appendix we evaluate the free energy and the internal energy for the two 
dimensional spin glass with a plus-minus measure. As stressed in Sec. (|12.2|) we get 
only a rough approximation, but it provides however a case in which an analytical 
calculation is possible. 

According to the general rule of Eqs. 1)22124(1 which, specifically, for a centered 
measure amount to Eqs. ((93194(1 . the free energy is readily obtained simply by using 
for ip Eqs. 1(102(1 and ((103() with the mapping 

ipi (tanh(/3J)) -^tp=^^i (tanh^(/3J)) . (B.3) 
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The internal energy U can be derived from the free energy F by using the formula 

Therefore, from Eq. H93() for the free energy per bond and per unit of couphng, 
e = E/{JNb), we have 

e = -tanh(/3J)-i^^, (B.5) 

where the factor 2 in the denominator of the second term takes into account that for 
any two sites we have a bond {Ni, = N/2). The derivative of ip can be calculated as 
follows 

dip dip dz 
'dp " 'd^W 



(B.6) 



dip 2z 1 p <9k , sin^ (uj) 1 

- ' K—duj — Ti (B-7) 



dz 1 + z2 TT /„ dz 1 I n ^2 o;v,2 ^ 2 ,.2o;„2 



1 + (1 - k2 sin^ (u;)) ^ (1 - k2 sin^ (uj)) ■ 



2 J , (B., 



9z 2 J tanh(/3 J) 
9^ ~ cosh^(/3J) 
whereas k is given by Eq. Ij 103(1 and its derivative is 

dn 1 — z^ 2^ 2 1 ~ m o^ 

(1 + Z2)2 ~ (1 + z2)2 ~ (1 + ^2)3 ■ ( 

By inserting the expressions (|B.6IB.8|) in Eq. (|B.5|) and by using for p the mapping 
l(B.3p . one gets the free internal energy e whose plot is reported in Fig. ((l()|l . 



Appendix C. Effective fields 

In our approach the key ingredient is the knowledge of the non trivial part of the free 
energy 0/ of the related Ising model, in terms of which the free energy is written as 

- PFi iVlog(2) + ^^i^^nlllog (cosh [PJ''^^)) + 0/ (tanh [PJ''^^)) ■ (C.IO) 

Once (pi is known, according to Eqs. ((23I24|) . in order to calculate </>, we have to 
implement in 0/ the following transformations 



and 

tanh 



jdfi tanh^ ipj,) ^ -fc)- + O (^^^ , (C.ll) 

(/3JW) Jdfi tanh iPJ,) = ^ + O (^^) . (C.12) 

As is well known, up to terms negligible in the thermodynamic limit, the fully 
connected Ising model with homogeneous coupling J^^\ has a mean field solution 
which, depending on the sign of j'^^ can be ferromagnetic or antiferromagnetic. It is 
convenient to distinguish the two sub-cases. 

Let us suppose that j'-^^ > 0. In this case the mean field solution for the related 
Ising model is given by 

-/3F/ =iVlog(2cosh(/3j(^)m/Arjj ^ (C.13) 
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where m/ satisfies the mean field equation 

mi = tanh (^fSmij'-'^N^ . (C.14) 
From Eqs. (|Cn?))l and (|C.13IC.14p . for the density ipi = (j)i/N wc then get 
If I (tanh {pj'-^^y) = log (cosh (/3j(^)m/7v)) - y log (cosh (/3J^^' 



(C.15) 

m/ =tanh(/3m/ J(-f > at) 

Equation IIC.15|) must now be used in the mapping of Eqs. 123I24|I . by using 
the transformations (|C.lllC.12p . To this aim it is convenient to use the relation 
cosh^(x) = 1/(1 — tanh^(a;)) and to rewrite Eq. IjC.lSp as 

1. . ox iV 



(tanh (Pj^^^y) = - - log (1 - m]) + — log (l - tanh^ (/3j( 
arctanh (tanh (/5J^^^)) m]N 



I) 



(C.16) 

m/— tanh (arctanh (tanh (/?J(^) ^ ^m / A^^ 

By inserting the transformations (|(/.11() and (|(/.12() in Eq. (|(].16p for N large we 
arrive, respectively, at 

fi i^j dp. tanh^ (/3J) j = -- log (1 - mlc) - , (C.17) 

and 

fj(^j d^, tanh(/3J)) =-ilog(l-r4)_^:^, (C.18) 

where the effective fields msc a-nd mp have been obtained, respectively, by using the 
transformations (jC.ll|) and (|C.12p and arc given by 



msG = tanh [[pjYmsG) , (C.19) 

and 

mp^ tanh {pJomp) ■ (C.20) 

The obvious interpretation of these fields is that, whereas m^? is the magnetization of 
the system (in presence of a disorder with variance J), mso represents ^qpA, ^-e., up 
to a square root, it is the Edward- Anderson order parameter, which below the critical 
temperature is non zero even for Jq ~ 0. 

Let us now suppose J^^' < 0. In this case the mean field solution for the related 
Ising model is given in terms of two effective fields ni'f^ and mf^ related to two 
sublattices a and h in which the given lattice A can be decomposed. The sublattices a 
and h are defined symmetrically so that the first neighbors of a site of the sublattice 
a are sites of the sublattice h and vice versa. In terms of the effective fields m^f' and 
mf^ the free energy of the related Ising model is given by 

-PFj = ^ log (2 cosh (/SJ^^'to^^^A^) 2 cosh (/JJ^^^mf ^ A^) ) 
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where rnf'^ and satisfy the mean field system 

TO^"^ - tanh f Pmf^ J^^^ N] , 

^ < (C 22) 

m^'^ = -tanh(/3m^")j(^)iVj . 

Similarly to what done in the previous case, from Eqs. (|C.10|) and ljC.21IC.22|l . for 

the density ipi we get 

ipi (tanh (pj^^^)) = \ log (cosh (/3 J^^^m^^^iv) cosh (/3 J^^^mf 'iv) ) 

+ ^^^"'f'"'' -^log (cosh (,.,,».)). ,C.23) 

and by inserting the transformations IjClljl and l|C.12p for N large we arrive, 
respectively, at 

1' 



ipi (^j dfi tanh^ (/3J)^ = - i log ((l - m^^^l^) (l - m 

(/3j)'m(f, 



) (fe) 

''SG'^^SG 



(C.24) 



and 



/ tanh = - 1 log ((l - m^^l') (l - m^^l^)) 



PJom^AF^AF 



(C.25) 



where the effective fields and m^gQ, and to^^ and m^p\ satisfy the following mean 
field systems, respectively 

- tanh ((/3J)2m('') 

\ (b) 
[ ^SG^~ 

and 



™SG = - tanh ((/3J)^m^"^ y , 



(C.26) 



m^AF ~ ^ tanh (pJofn'^AF I ' 
TO^]^. = — tanh (f3Jom,[^p ^ 



(C.27) 
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